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Abstract 

The tunneling rate, with exact prefactor, is calculated to first order in h for a closed 
FRW universe filled with perfect fluid violating the strong energy condition. The calcula- 
tions are performed by applying the dilute-instanton approximation on the corresponding 
Duru-Kleinert path integral. It is shown that a closed FRW universe filled with a perfect 
fluid with small violation of strong energy condition is more probable to tunnel than the 
same universe with large violation of strong energy condition. 
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1 Introduction 



The universe as a whole was the subject of intense investigations in the early 1960s with the 
work of DeWitt, Wheeler and Misner Lately, in early 1980s, the field underwent a rebirth 
to fully explain the very early history of universe. The idea of quantum tunneling, responsible 
for the birth of the universe, was first introduced by Atkatz and Pagels || and has been 
further developed by Vilenkin ||. The simplest possible quantum cosmological model is the 
spontaneous birth from "nothing" of an empty FRW universe with a constant vacuum energy 
density as its effective cosmological constant. Atkatz and Pagels showed that only a closed 
universe can arise via quantum tunneling. The rate of tunneling from "nothing" is obtained 
by WKB approximation for this model 0). 

The purpose of this paper is to investigate the issue of quantum tunneling for a closed FRW 
universe with a perfect fluid source violating the strong energy condition (SEC). It has been 
shown [|J that neither quantum nor classical wormholes, as instanton solutions "connecting two 
asymptotically flat Euclidean domains", exist for such a perfect fluid; however the violation 
of SEC is somehow reasonable in early universe || followed by processes such as quantum 
tunneling. So we are interested in finding the instantons as "tunneling solutions" for this type 
of perfect fluid. To the extent we are concerned with instanton calculations for tunneling rate, 
this model may effectively be replaced by an empty closed FRW universe with a cosmological 
constant. We shall calculate the tunneling rate by applying the dilute- instanton approximation 
to first order in h 0, on the corresponding Duru-Kleinert path integral 0. Its prefactor is 
calculated by the heat kernel method 0, using the shape invariance symmetry ||. 

In section 2, the Duru-Kleinert path integral formula and Duru-Kleinert equivalence of 
corresponding actions is briefly reviewed. In section 3, we introduce the cosmological model of 
closed FRW universe with a cosmological constant as a model without cosmological constant 
but effectively filled with a perfect fluid violating SEC. Finally in section 4, the tunneling rate 
is fully calculated to first order in h by applying the dilute-instanton approximation on the 
corresponding Duru-Kleinert path integral. 

2 Duru-Kleinert Path Integral 

In this section we briefly review the Duru-Kleinert path integral ||. The fundamental 
object of path integration is the time displacement amplitude or propagator of a system, 
(X b t b | X a t a ). For a system with a time independent Hamiltonian, the object (X b t b \ X a t a ) 
supplied by a path integral is the causal propagator 



Fourier transforming the causal propagator in the time variable, we obtain the fixed energy 
amplitude 



(X b t b | X a t a ) = 9(t a - t b ) < X b \ exp(-iH{t b - t a )/h)\X a > . 



(1) 




(2) 
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This amplitude contains as much information on the system as the propagator (Jff, t b \ X a t a ) , 
and its path integral form is as follows: 

(x b | x a ) E = J™ dt b J vx(ty AE ' h (3) 

with the action 

Ae = t b dt[^x 2 (t) - V(x(t)) + E] (4) 

Jt a Z 

where x denotes the derivatives with respect to t . In || it has been shown that fixed energy 
amplitude ([|) is equivalent to the following fixed energy amplitude, 

(X b I X a ) E = J o dS\f r (x b )Mx a ) j Vx(s)e iA ^ n ] (5) 

with the action 

A f E = f S ds{-^—x' 2 (s)-f(x(s))[V(x(s))-E]} (6) 
Jo 2f(x(s)) 

where f r and /; are arbitrary regulating functions such that / = fif r and x' denotes the 
derivatives with respect to time s. The actions A E and A E , both of which lead to the same 
fixed-energy amplitude (X b \ X a ) E are called Duru-Kleinert equivalent []. 

In the following section we shall use this equivalence to calculate the quantum tunneling 
rate. Thus we rewrite the action A E in a suitable form such that it describes a system with 
zero energy; as only in this sense can we describe a quantum cosmological model with zero 
energy Imposing E — in (^|), with a simple manipulation, gives 

A = | «/(.vwk^a-v) - vmm ( 7) 

where X denotes the derivative with respect to new parameter s' defined by 

s' = S'h (8) 
with 5* ClS db dimensionless scale parameter. 

After a Wick rotation s' = —it, we get the required Euclidean action and the path integral 

i> = /; ^(A-M)i 2W( ; (r)r v'H + v(a-m)}, ( 9) 

(X b I X a ) = J™ dSlMXMXJ I PA'(T)e-'»'"'] (10) 
where r is the Euclidean time. The action (H) is Duru-Kleinert equivalent of 

"6 , r M 



'0 



drf-X\r) + V{X{r))} (If 



where r a and r b correspond to t a and t b respectively, and X denotes the derivative with respect 
to Euclidean time r. 



1 Of course a third action A^ K Z is also Duru-Kleinert equivalent of Ae and A E but we do not consider it 
here ill. 
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3 Model 



We consider an empty closed FRW universe with a nonvanishing cosmological constant A. 
The system has only one collective coordinate, namely, the scale factor R. Using the usual 
Robert son- Walker metric we obtain the scalar curvature 



K = 6 



R 1 + R 2 
R + R? 



(12) 



Substituting (jTJ) into the Einstein-Hilbert action with cosmological constant leads to the 
action 

dt \A RR 2 + ^ R (i_ (13) 



with the constraint of Einstein equation 



4* 



0. 



(14) 



The system of ([13]) and (|TJ) is effectively equivalent to a model of closed FRW universe without 
the cosmological constant but with a stress-energy tensor of perfect fluid parametrized by 
constant vacuum energy density p vac and pressure p = —p vac [0 • 

Usually a perfect fluid is described by an equation of state such as p = (7 — l)p where the 
energy density p(R) and pressure p(R) may be functions of scale factor R. We then find, using 
the continuity equation, the behaviour of the energy density in a FRW universe M 



p(R)=p(Ro) 



(15) 



where R$ is the value of the scale factor at an arbitrary reference time. It is a well-known fact 
that SEC is satisfied only if 7 > |; however it is expected to be violated in the early universe 
H followed by such processes as quantum tunneling. So we are interested in investigating the 
perfect fluid for 7 < | and look for such processes as quantum tunneling for the corresponding 
cosmological model. Thus we rewrite the relation (|15D as 



p{R) = p{R ) 



Rq 



(16) 



where a > indicates the violation of SEC. As we shall see, to the extent we are concerned only 
with the action ( |I~3"D and the corresponding constraint (|14"D dealing with saddle point evaluation 
of the Euclidean path integral, we may substitute this energy density into the system (|13|) and 
(|I4T) as the effective cosmological constant A = p(R) 



dt 



-RR 2 + -R 
2 2 



R_ 

Rq 



(17) 



2 We take units such that |g 



1. 
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R 2 + 

where p(Ro) = ^7. Now, this system describes effectively a closed FRW cosmology in presence 
of a perfect fluid violating SEC. The issue of quantum tunneling for this effective model 
may be investigated in two ways: WKB approximation, and dilute-instanton approximation 
techniques. We shall follow the second approach here, in order to calculate the tunneling rate 

r. 



;is) 



4 Tunneling rate 



It is not suitable to apply the instanton calculation techniques to the Euclidean form of the 
action (|T7|). The reason is that the kinetic term is not in its standard quadratic form. It 
has been recently shown that in such cosmological model one may use the Duru-Kleinert 
equivalence to work with the standard form of the action. Using the same procedure, we find 
the Duru-Kleinert equivalent action in the cosmological model here as follows 



Tb 



dt 



(19) 



Now, the Euclidean action ( |T9"D has the right kinetic term to be used in instanton calculations. 
The Euclidean type Hamiltonian corresponding to the action ([T9|) is given by 



He = — —R 2 



Ro 



(20) 



whose vanishing constraint He = | gives a non-trivial instanton solution 

Ro 



R(r) 



(cosh(£))< 



(21) 



corresponding to the potential 



V(R) = -R 2 



1-1* 
Ro 



For R>0. 



(22) 



Each solution with a > (violating SEC) describes a particle rolling down from the top of a 
potential —V(R) at r — ► — 00 and R = 0, bouncing back at r = and R = Ro and finally 
reaching the top of the potential at r — > +00 and R = 0. 

The region of the barrier < R < Rq is classically forbidden for the zero energy particle, but 
quantum mechanically it can tunnel through it with a tunneling probability which is calculated 
using the instanton solution ([21]) . 



3 The constraint He = corresponds to Euclidean form of the Einstein equation. 
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The quantized FRW universe is mathematically equivalent to this particle, such that the par- 
ticle at R = and R = Rq represents "nothing" and "FRW" universes respectively. Therefore 
one can find the probability 

| < FRW{R ) | nothing > | 2 . 

The rate of tunneling T is calculated through the dilute instanton approximation to first order 
in h as || 

det'(-d 2 + V"(R)) 1/2 ^Mm MR) 1/2 
r = [ det(-d 2 + uj 2 ) ] 6 R [ ^h ] (23) 

where det' is the determinant without the zero eigenvalue, V"(R) is the second derivative of 
the potential at the instanton solution (pip, uj 2 = V"(R)\ji =0 with oj 2 = 1 for the potential 
(p2|), and Io(R) is the corresponding Euclidean action evaluated at the instanton solution (pl"|). 
The determinant in the numerator is defined as 

oo 

det>[-d 2 + V"(R)} = Y[\\ n \ (24) 

ra=l 

where X n are the non-zero eigenvalues of the operator — d 2 + V"(R). 
The explicit form of this operator is obtained as 

n- -2r d2 (« + !)(« + 2) , 2 , , . 

° = " h ^ crf^ + " ] (25) 



where we have used ( |21~| ) and (^) with a change of variable x = -. Now, in order to find 



exactly the eigenvalues and eigenfunctions of the operator (|25|) we assume a to be positive 
integer. By relabeling I = a + 1, the eigenvalue equation of the operator ( P5| ) can be written 

as 

A^,(x) = (£,-2Z + l)Vi(x) (26) 

with 

A,:=-^-^±^ + / 2 (27) 
fix 2 cosh 2 x 

where the factor a~ 2 is ignored for the moment. The equation ( ^6|) is a time independent 
Schrodinger equation. Now, by ignoring the constant shift of energy 21 — 1 and by introducing 
the following first order differential operators 



f Bi{x) := £ + /tanhx 
\ fi/(x) := -f + Ztanhx, 

the operator A/ can be factorized and using the shape invariance symmetry we have 



(28) 



Ms) = -^=B}{x)ik-i{x) (29) 
^_i(x) = -^ w B l (x)ii l (x) (30) 
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Therefore, for a given I, its first (bounded) excited state can be obtained from the ground state 
of I — 1. Consequently, the excited state m of a given I, that is ipi m , can be written as 



1pLm(x) 



2{2lU ' B}{x)BUx)...Bl +1 {x)-±— , (31) 



\ n™ =1 j(2/-j)2 m (m-l)! (V 7 m+iW cosh m x 
with eigenvalues E^ m = I 2 — m 2 . Also its continuous spectrum consists of 

B\{x) BU{x) B\{x) e** 



with eigenvalues E^ k = I 2 + k 2 where ipi k (x)ipi t k>(x)dx = S(k — k'). Now, we can calculate 
the ratio of the determinants as follows. First we explain very briefly how one can calculate 
the determinant of an operator by the heat kernel method 0. We introduce the generalized 
Riemann zeta function of the operator A by 

where A m are eigenvalues of the operator A, and the determinant of the operator A is given 
by 

det A = e-^ 0) . (34) 



It is obvious from the equations ([33|) and ([34]) that for an arbitrary constant c 

det{cA) = c u{0) detA. (35) 

On the other hand (a{s) is the Mellin transformation of the heat kernel G(x, y, r) Q which 
satisfies the following heat diffusion equation 

^,»,r) = -^, (36) 

with an initial condition G(x, y, 0) = 8{x — y). Note that G(x, y, r) can be written in terms 
of its spectrum 

G{x, y,r) = Y / e- x - T r m {^m{y). (37) 



An integral is written for the sum if the spectrum is continuous. From relations (p4[) and (|36D 
it is clear that 

I roc r+oo 

C A (s) = —- drr 8 - 1 dxG{x,x,r). (38) 

1 [S) JO J-oo 

Now, in order to calculate the ratio of the determinants in (|2"3"D, called a prefactor, we need to 
find the difference of the functions G(x, y, r) for two operators A/, Az(0), where 

AK0):=-^ + / 2 . (39) 



^Hcre r is a typical time parameter. 
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Considering the fact that A^ + 1 — 2/ (or A/(0) + 1 — 2/) has the same eigen-spaces as A; (or 
Aj(0)) and the eigen-spectrum is shifted by 1 — 21, we have 



G Al (o)+i-2i(x,y,T) = — e ^ (40) 

Z\/7TT 



l-l 



G Al+1 „ 2l (x,y,r) = £ ^ m (x)^, m (y)c- |(ro - 1)(a, - (B,+1))|T + 

m=0,m^l 



rOG 



dke-^r+^l^M. (41) 
In order to calculate the function ( Al +i-2h according to the relation fl38|) we have to take the 



trace of heat kernel G Al+ i- 2 i(x, y, t) where we need to integrate over \ipi tk \ 2 - Using the relation 
-%L=iJ) l k (x) = V»/-i,fc(a:) we have 



1 

dxif)l k (x)if>i >k (x) = - lim —= -J)* lk (x)^i_ x>k {x) + 

1 r+oo 

lim -77 — -^tki^i-iAx) + / ^i-iA^i-wi?)' ( 42 ) 



The first and the second terms appearing on the right hand side of the recursion relation ( f£2|) 
are proportional to the asymptotic value of the wave functions at oo and — oo, respectively, 
where the latter is calculated as 

-ikztm —ik±(m — l) — ik ± 1 exp(ikx) 



lim ibmk(x) = — 



71 



1 m —ol- 4- 9 



Substituting these asymptotic behaviours in the recursion relations between the norms of the 
wave functions i^ m ,k associated with the continuous spectrum (|32|) , then using the obtained 
recursion relations together with the orthonormality of discrete spectrum we get the following 
result for the difference of traces of heat kernels: 



dx(G Al+1 ^2i(x,x,r) - G Al ( 0)+1 - 2l (x,x,r)) = 

, ^/o, / ,u \ 1^ , f+°° exp(-((Z - l) 2 + k 2 )r 
£ exp(-|(m - 1)(2/ - (m + l))|r) - - £ m( / dfc PV VV 7 



(k 2 + m 2 )(k 2 + (l-l) 2 ) 
Hence, using the Mellin transformation (j38|) and the well-known Feynman integral 

1 



a(a 1 )T(a 2 )---T(a n ) J (t^ + t 2 D 2 + ■ ■ ■ + * n £) n )ai+« 3 +-+«» 1 ' 
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We finally get 



CA ; +l-2i( s ) - Ca ; (0)+1-22( s ) = 

E (\(m- 1)(2Z - (m + 1))|)- - ±1(1 + 1)(Z - iy^P(s + ~ ~) 

Z7T Z Z 



m=0,m=£l 



E m(Z - l)-( 2 *+ 3 )((/ - l) 2 - m 2 ) 2 F l( , + §, 1, a + 2, 1 - 



i r( s + |) ; „ 2(s+1) . 2 n , 3 , 1 _ o 1 (/-l 



,2 



^1x1-20^ + _ s + i, s + 2 , 1 - 



V^Fr(s + 2) 2 1V 2' Z 



(46) 



where /3 is the freta function. 
For s = 0, we obtain 

Ca,+i-m(s) - Ca ! (o)+i-2«(s)| s =o = -1. (47) 

This means that the operators A; + 1 — 21 and Aj(0) + 1 — 2/ have the same number of 
eigen-spaces (even though for both of them this number is infinite) since from the definition of 
Riemann's zeta function, it is obvious that its value at s = can be interpreted as the number 
of eigen spaces of the corresponding operator. The appearance of —1 on the right hand side of 
the relation is due to the ignorance of the eigen-functions associated with the zero eigen value 
of the operator A; + 1 — 21. Therefore, its number of eigen-states is the same as that of the 
operator A/(0) + 1—2/. In order to calculate the ratio of the determinant of the operators 
A; + 1 — 2/ and A/(0) + 1 — 2/ we need to know the derivative of their associated zeta functions 
at s = 0. Hence differentiating both sides of the relation fl4"7|) with respect to s and evaluating 
such integrals as 



L 



1 3 



dt(((l-l) 2 -m 2 )t + m 2 )-i logt = — — — — -(log2 + logm-21og(/-l+m)), (48) 

m((l — iy — wr) 



/ dt(((l - l) 2 - m 2 )t + m 2 r 2 log(((/ - l) 2 - m 2 )t + m 2 ) logt 
Jo 



ni 



(1 _ i _ m ) 

m ,\ 2 ^(logm-^- log(/ - 1) + 1 } )m - 2 log(/ - 1 + m)), (49) 

m\\l — \y — wr) / — 1 / — 1 

we get 

2(2/ - 1)! 

CA i+ i-2«(s) - Ca,(o)+i-m(s)I*=o = 1q s( (^_ 2 )!)2 - ( 50 ) 

Therefore, according to the relations fl23|), fl34"|), (|35|), fl4"7|) and (|50| ) the prefactor associated 
with the potential (p2"D , is calculated as 

1 (_d? _ 1(1+1) 



^((7^F(-S-^ + (^-l) 2 ))_((/-l)!) 2 



i , V , „ ^ " ~n7n~7 iTT - (51) 



det(j^(-£ + (l-l) 2 )) 2(2/ -1)! 
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Finally, the decay rate of metastable state of this potential is calculated as 



r = ^= J R 2 / " 1 exp(-^ ; °^ V2 ' ^ + 0(h) 



(52) 



where we have used the value of Jo at the instanton solution fl2~ID : 



/CO 
-co 




(53) 



The result (0) agrees with one obtained by the WKB approximation, and generalizes the 
decay rate obtained for I = 2 (or a = 1) in ||. 

Concluding Remarks 

In this paper we have calculated the tunneling rate, with exact prefactor, to first order in 
h from "nothing" to a closed FRW universe filled with perfect fluid violating the strong en- 
ergy condition. We were motivated by the belief that the strong energy condition is expected 
to be violated in the early universe, especially in such processes as quantum tunneling. The 
result obtained here somehow confirms this fact in that the tunneling rate ( |52"D increases for 
higher values of the positive integer a (or I) such that for a — > oo the very high tunneling rate 
corresponds to the point 7 = 2/3 below which the strong energy condition is violated 0. On 
the other hand for a = 1 the tunneling rate has its minimum value corresponding to 7 = 
(cosmological constant |4[]) as more violation of strong energy condition. In this respect, the 
main result is that a closed FRW universe filled with a perfect fluid with small violation of 
strong energy condition is more probable to tunnel than the same universe with large violation 
of strong energy condition. 
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